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It is known that in certain non-bipartite quasi-one dimensional spin systems in a magnetic 
field, in addition to the usual Pauli coupling of the spins to the field, new parity breaking three 
spin interactions, i.e. chiral spin interactions, are induced at higher order due to virtual processes 
involving the intrinsic electronic nature of the underlying spins. The strength of these interactions 
depend strongly on the orientation of the field, a feature which can be exploited to detect chiral 
effects experimentally. In many spin systems, these chiral interactions are generated and should be 
taken into account before any comparison with experiments can be made. We study the effect of the 
chiral interactions on certain quasi-one dimensional gapped spin half systems and show that they can 
potentially alter the physics expected from the Pauli coupling alone. In particular, we demonstrate 
that these terms alter the universality class of the C-IC transition in spin-tubes. More interestingly, 
in weakly coupled XX zig-zag ladders, we find that the field induced chiral term can close the singlet 
gap and drive a second order transition in the non-magnetic singlet sector, which then manifests 
itself as a two component Luttinger liquid-like behaviour in the spin correlation functions. Finally, 
we discuss the relevance of our results to experiments. 



I. INTRODUCTION 



Quasi-one dimensional quantum spin systems have long been known to exhibit a spectrum of behaviours, ranging 
from the gapless critical behaviour in spin-i integer systems, to gapped spin Peierls systems and Haldane phases in 
integer spin systems. The effect of an external magnetic field is usually described through the Pauli coupling (h • S) of 
the spins to the field. These fields can have dramatic effects on the spin systems. For example, in gapped spin systems, 
the magnetic field induces the commensurate-incommensurate (C-IC) transition, where the IC phase is characterized 
by gapless spin excitations described by a Luttinger liquid and the corresponding spin correlation functions exhibit 
power law behaviours whose exponents, however, vary from system to system. However, it is important to remember 
that in real materials, the spins are electrons whose charge degrees of freedom are frozen by strong electron-electron 
interactions. In a physical system, the magnetic field couples not only to the spin degrees of freedom, via the Pauli 
term, but also to the orbital degrees of freedom of the electrons. This results in non-trivial additions to the effective 
spin hamiltonian. In fact, it is known that in the limit of large^Coulomb repulsion, this coupling can in principle 
induce a three spin chiral interaction between neighbouring spinsEln. Such a term arises in two dimensional systems on 
the triangular and Kagome lattices. These chiral terms also arise in certain quasi one dimensional systems like zig-zag 
ladders, frustrated ladders, spin tubes, saw tooth chains, etc, whose geometries permit the formation of a closed 
plaquette involving three neighbouring spins. The possible generation of such chiral terms in quasi-one dimensional 
spin systems stemming from the underlying electronic nature of realistic models has been systematically overlooked 
in the literature. In this paper, we show that these chiral terms are indeed relevant and can change the behaviour of 
certain spin systems like spin-tubes and zig-zag ladders in magnetic fields. The possible effects engendered by these 
terms should be studied before sensible comparisons with experiments in a magnetic field can be made. Another 
important motivation is that the effects of these chiral interactions depend strongly on the direction of the applied 
field, and could be detected experimentally in measurements of various quantities ranging from specific heat, NMR 
rates and dynamical structure factor measurements. 

One class of systems, where the chiral term will surely be relevantpace systems which support chiral excitations 
in zero field. Examples are the recently studied zig-zag spin s ladderset j where certain parity breaking interactions 
induce spin nematic ground states in zero magnetic field and spin tubeaj, whose ground states have chiral properties 
in the absence of a magnetic field. These spin tubes were also shown to exhibit an interesting C-IC transition where 
the IC phase is described by a two component Luttinger liquid describing gapless spin and chiral excitations. In these 
systems, the new field induced chiral term can completely alter the physics expected from the Pauli coupling alone. 
Normally, since we do not know of physical probes which couple directly to the chirality, it is quite hard to perceive 
experimentally the chiral aspects of the ground state. However, the competition between the chiral term and the 
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Pauli coupling of the spins to an external magnetic field can be exploited to obtain experimental signals of chirality 
in certain systems. The main objective of this paper is to study the influence of these chiral terms on the gapless 
incommensurate phases of certain spiflpSjystems. These are especially pertinent given the surfeit of experiments on 
spin systems in strong magnetic ficldstllij. 

Thcj_paper is organised as follows: we first outline the derivation of the chiral spin term starting from the Hubbard 
modetl. We then study the effect of this term on the spin-tube and other quasi-one dimensional systems, in the limit of 
strong planar interchain couplings. We then discuss the very interesting example of the effects of the chiral interaction 
on the weakly coupled zig-zag ladder, where the magnetic field induced chiral interaction triggers a transition in the 
spin singlet sector. Finally, we conclude with a discussion of the experimental consequences of this interaction. 



II. CHIRAL INTERACTION 



In the presence of an external magnetic field, it is customary to only consider the Pauli interaction. However, this 
completely ignores the coupling of the field to the orbital part of the motion which can result in non-trivial corrections 
to the effective spin hamiltonian. We start by considering a system of interacting electrons described by a Hubbard 
model 

H = - Uj [c\cj + h.c] + n^mi , (1) 

(ij) * 

where the first term describes the hopping of electrons with spin a from sites i to j and the second, the on-site Coulomb 
repulsion between electrons. The Pauli coupling of the electronic spins S to the external magnetic field h is given by 
/ioffh • S, where fj,o = \e\h/2mc is the Bohr magneton and the gyromagnetic ratio g ~ 2. Hereafter we set g/io = 1. 
The magnetic field also couples to the electronic motion via the hopping matrix elements Uj , which become complex, 

i.e., it picks up a phase Uj —> Uj exp[ie f. A • dr], where, A is the vector potential associated with the external field h. 
If the electron hops around a closed loop, we immediately see that the net hopping matrix element is proportional to 
the magnetic flux enclosed by the loop. Since an external magnetic field explicitly breaks time reversal invariance, it is 
interesting to study whether time reversal violating terms are generated in the spin hamiltonian in the large U limit. 
At second order in perturbation theory, no parity breaking terms are generated and one obtains only the Heisenberg 
term Si • Sj. Chiral terms, can therefore, be obtained only in higher order perturbation theory. An example of a chiral 
term which respects all spin symmetries and yet violates time reversal and parity, is the triple product of three spins 
Si.Sj x St. Is such a term generated within the Hubbard model with nearest neighbour hopping t in the presence 
of an external magnetic field ? At half filling, it can be shown that no parity and time-reversal violating terms are 
induced on bipartite lattices!! However, on non-bipartite lattices, such a term is explicitly generated at order t/~ 2 in 
perturbation. 

24^3 

H ch irai = -jp sin(e$) ^ S i- S J x S * ( 2 ) 

(ijk) 

where the summation is over sets of nearest neighbour sites which form triangles. We have restricted ourselves to 
nearest neighbour hopping t in the above. $ is the flux enclosed by the plaquette formed by the sites i,j and k. In 
terms of the magnitude of the field h, $ = Ah cos(#) where A is the area of the plaquette and 8 is the angle made by 
the magnetic field with respect to the normal to the plane of the plaquette. Substituting this in (|^), we see that no 
chiral term is generated for magnetic fields applied parallel to the plane of the plaquette. Additionally, the strength 
of the chiral term varies from a minimum to a maximum as a function of the direction of the applied field. For real 
materials, the typical area enclosed by these triangles and hence the enclosed flux being very small, the sin(e$) term 
in (J2j) can be expanded to yield 

24^3 

Hchirai = -jj^eAhcosO Si.Sj x S fe (3) 

(ijk) 

Since the chiral term is generated at higher orxler, its strength is lesser than the Heisenberg exchange as well as 
the Pauli interaction. Despite this, it was showncl that in certain cases the chiral term may have a larger effect on 
the ground state than the Pauli coupling of electron spins to the external magnetic field, thereby drastically affecting 
the ground state properties of the model. One example is the case where the ground state of the system in zero field 
is doubly degenerate with opposite chiralities. Then, when an external magnetic field is applied, in addition to the 
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effects of the Pauli term, the chiral degeneracy is broken by (||). The effect could be more pronounced in the case 
where the ground states are degenerate chiral singlets because here, though the chiral term still breaks the degeneracy, 
the Pauli interaction term has no effect on the ground state for fields smaller than the singlet-triplet gap. 

Since, no chiral interactions arc generated for fields applied parallel to the plane of the triangular plaquettes, 
the behaviour of the system could, in principle, be different depending on whether the field is applied parallel or 
perpendicular to the plane of the system. This is especially relevant for the critical behaviours in the IC phase, where 
exponents could show a variation with the field direction. Such a behaviour, provides us with a physical handle for 
determining whether an isotropic spin system is quasi-one dimensional or not i.e., if the behaviour is different for 
the two directions of the field, one can presume that this difference stems from chiral interactions which necessarily 
require the system to be atleast quasi-one dimensional! 



III. STRONGLY COUPLED SPIN TUBES AND LADDERS 



In this section, we discuss the effects of the chiral term on the IC phase of certain strongly coupled spin-gap systems, 
whose geometry allows for the existence of chiral terms in a magnetic field. In the absence of chiral inter^tkma, the 
behaviour of these systems in the IC phase, has been well studied and catalogued in the literatureBHO As 
we will show below, the chiral term does not necessarily alter the nature of the incommensurate phase always. The 
effects of the chiral interaction are drastic or trivial depending on the nature of the ground and excited states of the 
system. 



A. Spin tubes 

The first example we consider is the strongly coupled spin tube (i.e. a three-leg ladder with periodic boundary 
conditions) studied in Ref. ||. As already mentioned spin tubes, exhibit a new class of C-IC behaviour. Starting 
from the limit of strong-rung coupling (see Figj^), the low-energy effective Hamiltonian (LEH) describing the gapless 
incommensurate phase with finite magnetization was derived in Ref. ^. This low energy Hamiltonian was that of 
a two-component Luttinger liquid describing gapless spin and chiral excitations in the system. The existence of a 
two-component Luttinger liquid phase has important physical consequences. Namely, it implies that both chiral and 
magnetic excitations contribute to a non-zero magnetic susceptibility and a strongly enhanced T-linear specific heat. 
This different behaviour stems basically from the presence of a degenerate ground state with different chiralities, within 
an isolated triangle. This immediately poses the interesting question, what is the effect of the chiral interaction (|^) 
on the two component Luttinger liquid behaviour? Clearly, we expect the chiral term (^) to induce dramatic changes, 
because it will break the degeneracy of the ground state. To see this, we derive the low-energy effective Hamiltonian 
for the spin tube in the presence of the chiral interaction (J3J) . We will use this effective theory to discuss the physics 
of the IC phase and the resulting spin-spin correlation functions. 

The Hamiltonian of the spin tube in the presence of an external magnetic field h and a spin-chiral interaction is 
given by, 

N 3 N 

H = }^ ^(JS; iP • S l+1 . p + J ± ■ S lyP ■ S. t . p+1 -h ■ Si, p ) + /i c h^ S;,i ■ (S ii2 x S ij3 ), (4) 

i— 1 p—1 i—1 

where (i,p) denote the site and chain indices respectively and periodic boundary conditions imply that the site (z, 4) 
is identified with the site (i, 1). The first two terms represent the spin Hamiltonian for the spin tube, the third the 
Pauli coupling of the spins to the magnetic field and the last term is the chiral term derived earlier. J is the coupling 
along the chain, J± the transverse coupling between the chains and, \i c = ^-eA cos 9 is the strength of the spin-chiral 
interaction. As in Ref. ||, we consider the limit J± J. In this limit, the properties of the system can be studied 
by perturbing in J around the limit of decoupled triangles. For J — 0, the system consists of independent triangles. 
The eight states of a given triangle split into two spin-^ doublet ground states, with energy E s — — 3Jj_/4 and an 
excited spin-| quadruplet, with energy E t — 3Jj_/4. The spin doublet ground state has chiral eigenvalues ±. For 
strong magnetic fields applied parallel to the plane of the triangle, no chiral terms are generated. Consequently, the 
degeneracy in the ground state is only partially lifted and one reverts back to the case studied in Ref. ||. On the 
other hand, for fields perpendicular to the plane of the triangle, the ensuing Pauli and the spin-chiral terms lift all 
the degeneracies. For instance, for the spin doublet states, the Pauli and chiral terms contribute ±h/2 and ±/x c -\/3/4 
respectively, to the energy . Note, however, that the chiral term has no effect on the ferromagnetic spin | excited 
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states. This is expected since the chiral term, which is nothing but a triple product, has a non-zero expectation value 
only for non-collinear spin configurations. Therefore, for the field applied in the z direction, the resulting energies are 



E t = E t - hm 3/2 

E s = E s - hmt/2 T — | — 

Here m s just refers to the S z value of the state considered. There exists a critical field h c , such that for h < h c 
the ground state has magnetisation m = i and for h > h c the ground state has m = |. In the absence of chiral 
interactions, the critical field h c required to close the gap is given by h c — E t — E a = 3 Jj_/2. However, in the presence 
of chiral interactions, the critical field is slightly increased and is given 

= (E t - E s ) _ 3Jj_ , , 

1 4 Z 2 

We therefore see that the chiral interaction changes the magnitude of h c in a subtle way. Note also that the magnitude 
of the critical field depends on the direction of the applied field. It is minimum for a field applied in the plane of the 
triangle i.e., x direction where h c = 3Jj_/2 and is maximum for the field in the z direction. A small non-zero value 
of the coupling J, broadens the abrupt transition from a state with m = ~ to m = |. This results in two critical 
fields H c i (where the m = h state stops being the ground state ) and H C 2 (where the magnetisation saturates at 
m = § corresponding to a ferromagnetic state) with the magnetisation varying smoothly in between. The difference 
H C 2 — Hex is of the order of the coupling J. We now write an effective hamiltonian describing the low energy physics 
of the spin tube close to the critical magnetic field. Since at h c the m s = § state becomes degenerate with only one 
of the chiral states with m s = i we retain only these two states per site : 

as = | TTT) (6) 

«i = ^[im>+iiuT)+i a im>] (7) 

where j = exp(27ri/3). We now introduce a pseudo-spin half operator T which acts on the subspace containing the 
above two states. By inspection, we can rewrite the original spin operators S in terms of T i.e., 

= ( 8 ) 



Sl P = -(l+T?) (9) 



Using this, the effective hamiltonian, to first order in the coupling between the triangles is now given by 

N N N 

Heff. = g £( T i +T <+i + T-T+ +1 ) + - £ T?Tfa - h eff . £ Tf (10) 

i—1 i—1 i—1 

where h e f / = h — h c — ^ and contains the effects of chirality. The previous derivation is equivalent to that of the 
low-energy effective Hamiltonian for three coplanar spin- ^ chains considered in Ref . [TJ|. The model ( |l0| ) is the simple 
XXZ chain in a magnetic field. Since the coefficient of the z component of the interaction is smaller than that of the 
transverse components, we see that the above system corresponds to that of the XXZ chain in its gapless phase. A lot 
is known about the model and results for the exponents of various correlation functions as a function of the applied 
field can be obtained from Ref. [?]. 

Using the Jordan- Wigner transformation we fermionize ( |To| ) , to obtain 

#eff. = ^ ^2i c i c i+i + h.c] + mrij+i ~ Li^2ni + const. (11) 

i i i 

H = h — h c — J/3. The values of the critical fields H c \ and H C 2 for the spin-tube can now be obtained from the 
effective model (pd|). When h = H c i, the band of spinless fermions is empty. This occurs when the chemical potential 
\x equals the band minimum i.e., when /i = —J which in turn yields H c \ — h c — 2J/3. Similarly, H C 2 corresponds 
to a completely filled band of electrons or equivalently an empty band of holes. The simplest way to describe this 
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situation is to perform a particle- hole transformation on the Hamiltonian (|TT|): c\ — > hi. Up to a constant, the new 
Hamiltonian reads: 

#eV = ~ Y^Hhi+i + h.c] + ^ YJ n * n^i "WE ( 12 ) 

where the chemical potential Hh = —I 1 + 2 J/3 . In terms of holes, H C 2 corresponds again to the chemical potential 
where the band stars to fill up, thus we find H C 2 = h c + 2 J/3. 



The Hamiltonian (10) can be easily bosonized to obtain the various correlation functions and exponents in the 



incommensurate phase. Using the standard expressions for bosonization 

T z {x) = o x <p+ (-1)« (13) 

7r 7ra 

T+(a;) - - 7 ^e~ 4e [(-l)» +cos20(x)] (14) 
\Z2na 

we obtain the following continuum theory 

U = h I dx ^ uK ^ n ) 2 + (^)( 9 ^) 2 ] + \j dxh eff d x ^ (15) 



Hp 1 — j - 

where n is the momentum conjugate to the field cj) where u = % _ u jf\ is the spin wave velocity, and 



2K=[l--cos- 1 ('^)]- 1 (16) 

This indicates clearly that in the presence of chiral interactions, the spin tube in the IC phase, behaves as a one- 
component Luttinger liquid, describing spinless excitations alone. In the absence of the chiral term, both spin and 
chiral excitations were gapless which then led to the two component Luttinger liquid behaviour. Here, we see that 
the effect of the chiral term is to open a gap for the chiral excitations, which then do not contribute to the low energy 
behaviour. The magnetic field term, described by the gradient in (|l5|), can be eliminated by a simple shift of the 
(j> field i.e., (j> — * <f> + ■nmx where the magnetisation m oc h e ffK/u. For small values of the field /i e //, the effective 
magnetisation m increases linearly with h. This field shift results in the appearance of incommensurate modes in the 
spin-spin correlation functions. Note that K varies from K = 1 for the case of the XX- antiferromagnet to K = | 
for the isotropic XXX system. For the present problem, J z = J/3 resulting in K = 0.83 and u = 1.2J. It is well 
known that for this value of K the system is gapless and lies in the XXZ universality class. However, due to presence 
of marginal operators in the theory, the Luttinger parameter K now varies with the magnetisation. This variation 
depends also on the anisotropy J z l J = 1/3 and can be obtained numerically using the Bethe ansatz equations of Rcf. 

1 . 

In the critical incommensurate (IC) region characterized by a finite magnetization, we can easily calculate the 
various correlation functions of the effective spin chain and hence that of the original spin tube 

x ; z = (s;( x ,t)s;(o,o)) (17) 

X± = (S+(x,t)S- (0,0)), (18) 

where p is the chain index and ± refers to the intra and interchain correlation, respectively. From (^|), we see that 
the magnetisation m in the spin tube is related to the magnetisation of the effective system by rh = (1 — m)/3. The 
first correlation function is useful for neutron scattering experiments, whereas the correlation function ( p~8| ) defines the 
staggered susceptibility that is useful for the calculation of NMR relaxation rates. Using (||) which relate the original 
spins to the pseudo-spin i variables, we obtain 

X zz - \ [(T z (a;,i)T*(0,0)) + (1 + (T z )) 2 ] (19) 
X+- = I (T + (x,t)T- (0,0)) (20) 
which in turn leads to the following result 
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(1 - m) + - cos7nr(l - 2m)(x - ut) + const K. , , , . 

y J 9 y Jy ' \(x-ut) 2 (x + ut) 2 J 

(_)*/a( x 2 _ u 2fy^ + const coa (2Trmx)(x 2 - u 2 t 2 )-^ +K -^ 



(x - ut) 2 (x + uty 



(21) 



(22) 



where, u is the spin mode velocity. The correlation functions so found show a power-law decay in space and time. As 
anticipated, we see that two incommensurate modes appear at the wave vectors Q = 27r(l — m), for \ zz an d Q = 2wm 
for Xp~ i resulting in a different behaviour in the plane perpendicular to the field and along the field. The exponents 
of these power-laws now vary with the magnetisation m. These can be obtained by numerically solving the Bethe 
ansatz equations of Ref. |[ As the magnetisation m — > 0.5 or equivalently when the average number of fermions goes 
to 1, the fermion band is completely filled and we expect to recover the physics or free fermions whose K = 1. 

The behaviour found for the spin-spin correlation functions in the present case, is wholly different from that found 
in Ref. | 



(u a ty 



cos(Qi(x, t)) + cos(Q 2 (x, t)) 



cos(Q 3 (x,t)) 



{Uatf 



{u b tf 



x 2 - (u a t) 2 
+ cos(Q 4 (izi, t)) 



(u b ty 



(uaty 



(u b tf 



(23) 



Here, u a , Ub are the velocities of the spin and chiral excitations and K a , Kb are the Luttinger parameters associated with 
these excitations. Qi. Q2, Q3 and Q4 are functions which fix the incommensuration. One of the primary differences 
between (^3|) and (^2|) is that the former shows an explicit dependence on two different velocities u a and it& related 
to the magnetic and chirality modes and exhibits incommensuration at many more wavevectors. These velocities 
were different and in addition to the sum of the power law correlations exhibited by the two components, there is an 
interference between the two which contributes a third term with a different exponent in the correlation function for 
the spins. The two components will be distinguished when an extra field which couples to the spin current in the 
transverse direction is applied. 

To summarize, when the field is applied parallel to the plane of the triangles, no chiral interaction is induced and 
the system exhibits the two component Luttinger behaviour of (p3|). However, when the field is applied along the 
length of the tube i.e., perpendicular to the plane of the triangles, the chiral interaction plays a dominant role and 
we recover the single component behaviour indicated in ( |2l] , p2[ ). The exponents of the various power laws in the 
correlation functions will also be different in the two cases. Moreover, the value of the critical fields H c \ and H C 2 are 
also different for the different field directions. These predictions are strong signatures of the chiral term and they 
could be verified experimentally. Moreover, in the past there have been a lot of discussions of how to experimentally 
establish the existence of chirality in systems. In spin tubes and other systems, the change in critical behaviour can 
be perceived as an indirect test of such chiral exchanges in realistic systems. 



1. Frustrated spin ladders 



We hpee also analysed the effect of these chiral interactions on certain strongly coupled frustrated ladders considered 
by MilaO and shown in FigJ|. In all of these cases, we find that the chiral term has no effect whatsoever on the C-IC 
physics seen in these systems. This is due to the fact that in each of the examples considered, the unperturbed ground 
state (equivalent of the spin ^ states on the triangles in the spin-tube) has no interesting chiral properties i.e., they 
have no low lying chiral modes. As shown above, we introduce pseudo-spin operators to derive the relevant effective 
hamiltonians close to h c . In all these models, it turns out that the chiral term is effectively zero when written in terms 
of these operators. Consequently, the physics of these strongly coupled frustrated ladders is completely unaffected by 
the chiral term. 



IV. WEAKLY COUPLED ZIG-ZAG LADDERS 



In the previous section, we considered the effect of the chiral interaction on the spin-tube as an example of a strongly 
coupled quasi-one dimensional systems. Here, the possible non-trivial effects generated by the chiral term could be 
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deduced by merely studying the the chiral properties of the relevant "unperturbed" ground state in the absence of the 
external field. Such an approach fails for weakly coupled spin systems like the zig-zag chains and frustrated ladders. 
There, to begin with, one lacks a schematic picture, of the ground state(s) which is highly correlated. Neverthless, 
in such cases, one can use bosonization methodsE30 to study the effects of the chiral terms on the long wavelength 
physics of these systems. Here, we consider the example of the zig-zag ladder, where the chiral terms are relatively 
simpler than those generated in frustrated ladders. We will restrict ourselves to the limit of very strong anisotropy 
that corresponds to two coupled XX chains. (We assume that the anisotropies are generated by spin-orbit couplings, 
which do not affect the form of the chiral interaction.) The XX zig-zag ladder was studied in Ref. || using a mean-field 
treatment of the bosonized hamiltonian. Here, we analyze the effect of the chiral term on the mean-field solution via 
a similar treatment. 

The zig-zag ladder ( Fig.^j) has spins S on .chain 1 and spins T on chain 2. The bosonized hamiltonian for the zig zag 
ladder has been derived by Nersesyan et alJ3. The bosonic fields 0i and 02 and their duals 6\ and 62 are introduced 
to describe the spins in the two chains, S and T. In terms of the the symmetric (antisymmetric) combinations, 
0s, a — (01 ± 02)/v / 2 _ and 9 St a = {0\ ± 6 ) 2)/v / 2 _ which describe the triplet (singlet) sector, the hamiltonian describing 
the long wavelength physics of the anisotropic XX zig-zag ladder is given by 

H = H + Hi n t + H mag (24) 

where, 

a—s, a 

H in t = -yd x 6 s sin V26 a (25) 
and the Pauli coupling of the spins to an external field h (which affects only the triplet sector) is given by 

H mag = fl^-d x (f) s (26) 
7T 

The spin-wave velocity v cx J a and 7 is the strength of the TS%-7,&!g coupling. Note that the interaction term in J25"| ) is 
marginal. Due to the anisotropy, we only consider the case where the field is applied along the z-axis perpendicular 
to the plane of the zig-zag ladder. A magnetic field applied in the plane of the zig zag chain is very complicated, since 
even the simple Pauli coupling generates a gap in the triplet sector. In the absence of any field, the ground state 
of this model was shownQ to be a spin nematic with gapless triplet excitations and gapped singlet excitations. This 
state was also characterized by a power law decay of the staggered transverse spin correlations at an incommensurate 
wave vector Q — tt — S, where S depends on v and 7. A Pauli coupling to the magnetic field in the z direction has no 
effect on the staggered part of the correlation and induces only a trivial incommensuration in the uniform part of the 



spin correlation function. To the hamiltonian ( |24| ) we now add the chiral interaction term (27) involving the spins on 
every triangular plaquette. It takes the form 

(27) 

Hchiral = Hch{Si — Tj+l) ■ Tj X Sj+i (28) 

for every pair of spins Si, Ti. In terms of the continuum fields, 

-^sinV20 o (x)(l + 2v^a0 B o )+ (29) 



Hchiral = —(J, c h I dx 

1 



{2-k} 



;{d x 4)ad x 9 s - d x cj) s d x 8a) 



(30) 



We see that the first term in (|29|) is a relevant operator of dimension one that can have the effect of eliminating the 
gap in the singlet sector. The other terms are all marginal. The first two terms act exclusively on the singlet sector, 
whereas the last two terms couple the singlet and triplet sectors. Moreover, these terms will compete with H int in 
(^||). On physical grounds we speculate that the coupling of the chiral interaction to the singlet sector arises from 
the hidden chirality in the ground state. The presence of so many competing operators renders the problem highly 
complex. Here, since we are interested only in the qualitative effects of the chiral term, we restrict ourselves to a 
simple self-consistent mean-field analysis which leads us to the following Hamiltonian 

Hmf = H s + H a (31) 
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with 

h s = v - [(d x e s f + (d x <p s ) 2 } + (h + (3Lz)d x ct> 8 + {jKt - (3K 2 )d x e s , (32) 

H a = V - [(d x 9 a ) 2 + (d x (b a ) 2 ] + (ai + a 2 K 2 + jL 2 ) sin V29 a + 

+ (a 2 K 1 - 0L 2 )d m <f> a + /3L±d x 9 a , (33) 

K 1 = {wa.y/2,e a ) K 2 = {d x <f> a ) (34) 

Li = K 2 = {d x <t> a ) L 2 = (d x 9 s ) (35) 
L 3 = (d x 6 a ) U = (d x s ) , (36) 

and a.\ = —^ c h/2i: 2 , a 2 = — \/2~n c h/ir 2 ,/3 — /j c /2tt 2 . The triplet channel (s) can be solved easily by eliminating 
the terms linear in 9 S and <j) s through the field shifts: 9 S — > S + (7^1 — /3K 2 )x/2, and — ► <j) s + (ft + (iLz)x/2. 
Analogously in the a channel, the gradient of <j) a can be eliminated by a simple shift 4> a — > a + (a^-fta — f3L 2 )x/2. 
Self-consistency of the solution then leads to the following relations 



where the mean-field parameters 



K - = 2 ^2-4 7 Xl ' L2 = 5 (/3A - " " /V 1 ' ' '" ,7 1 



L 3 = -i(ft + / 9L 3 ) (38) 

Here we have assumed that the ground state of the system is found in the sector with nonzero spin current d x 9 s . 
^From the very structure of the mean field hamiltonian (|33"|), we see that though triplet excitations remain gapless, 
the situation is not so clear for singlet excitations. The hamiltonian in the singlet sector takes the form 

H a = \ [{d x B a f + (O^a) 2 ] +^smV26 a + (3L 4 d x 6 a (39) 

where the coefficient of the sine term which generates a gap is given by 

* = ai + ^4 Kl ' ( 4 °) 

Note that the structure of H a is similar to that of systems with gapped triplet excitations in an external magnetic 
field ( the equivalent of the "external" magnetic field in this case being the mean field parameter PL4). The bare 
gap in this case is A oc /1 2 / 3 . ^From the structure of the mean field equations, we see that since fi depends on the 
external magnetic field, the gap is a varying function of h. The field L4 tends to reduce this gap. Moreover, it is 
known that these systems exhibit the C-IC transition when the strength of the field attains a critical value which 
equals the magnitude of the gap. For fields greater than the critical field h c (not estimated here), the system has 
gapless excitations at incommensurate values of the wave vector. We, therefore, expect within mean field two kinds of 
solutions: in the first one the gap induced by the sine term survives and, consequently, there is no net "magnetisation" 
i.e., (d x 9 a ) = 0. This implies that K\ 7^ and L3 = 0. Using (|38|), we see that L4 = —h/2 in the gapped phase. 
The second solution describes the "incommensurate" phase with a finite "magnetisation" m a i.e., L3 = (d x 6 a ) — m a . 
In this phase, the effective field in the singlet sector is given by L4 — — + f3m a ). One or the other solution will 
be energetically favoured as the field is tuned. In the gapless IC phase, the gradient in 9 a can be absorbed into the 
quadratic part by the shift 9 a — > 9 a + (iL^xj2\ the singlet sector is now effectively a Luttinger liquid described by 



ti a = — - — 



K a {h){d x( j> a f + —^—{dJaY 
K a (h) 



(41) 



where the new velocity and the Luttinger liquid parameters v a ,K a both vary with the field. The results obtained 
earlier in the context of the normal C-IC transitiorO'Efl can be used to obtain K a {h). In fact, K a (h = h c ) = 2. 
In conjunction with the gapless spin excitations described by H Sl this results in a two component Luttinger liquid 
like behaviour! The possibility of such a transition, however, depends crucially on the mean field solutions. The 
calculation of the regime of validity of the two solutions requires the knowledge of the exact expectation value of the 
mass term (siny^^o) as a function of the prefactor of the "magnetic field" like term d x 9 a - This expectation value is 
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rather difficult to estimate and has only been calculated either for zero magnetic field orjn-the limit of fields larger 
than the gap opened by the sine-Gordon term i.e., in the gapless phase where A//3L4 <C lll3'E2l. Here, we assume that 
the chiral term has indeed driven the singlet sector across the transition into the gapless "incommensurate" phase and 
solve for the mean field equations. In this limit, 

sinV26» a \ = cf3L 4 n (42) 



where c is some constant and /i given by (40 ) is a function of the other mean field parameters. Using ([42| ) and 



we can now solve explicitly for K\ to obtain 



01 



Kl f-^-^M + ^h + i 2 ) (43) 

fl = U\ + 



/-i(M(/3 2 - 4) - 2( a 2 + a2/ 3 7 + 7 2) • 

where f(h) = (c/i) -1 . We find that the mean field solutions are self-consistent with the limit A/PL4 <C 1 used to 
obtain them. This limit translates to the condition h 3> and allows for a finite h c . 

In conclusion, we have shown that, while the Pauli interaction with the field leaves the physics of the XX ladder 
completely unchanged, at the mean field level the chiral term radically alters the physics of the zig-zag chain. It 
generates a transition to a gapless regime in the singlet sector for large enough values of the magnetic field. The 
gapless singlet excitations manifest themselves in specific heat measurements and more importantly, in the spin 
correlation functions. A new peak appears in the specific heat as a function of the magnetic field. Such a transition 
in the singlet space will result in a change of the exponents of power laws in the transverse spin correlation functions, 
as we will see below. 



A. The correlation functions and relevance to experimental quantities 

We conclude the analysis of the effect of the chiral interactions on the zig-zag chain with the calculation of the 
spin-spin correlation functions. They could be useful for the neutron scattering experiments and the NMR relaxation 
rates measurements. ^From the previous mean-field solution we get the following shifts of the fields 

9.^9. + *^* 
2v 

9a^9a + f-X 

2v 

<Pa^ <Pa + -^~ X 

j. 5 ' 

<Ps —> <Ps + TT X 

2v 

where 

7 (/3 2 - 4) - 0(a 2 + g/27) 
P= —± 

= a 2 + — 

for all values of the field h and 

p' = 0; 5' = h for h < h c 

p' = — — (h + /3m a ); 8 1 = h + 0m a for h > h c 

In the C phase (h < h c ), with gapless triplet excitations and gapped singlet excitations, (8 S ) = and 9 a has a non-zero 
expectation value determined by the position of the minima of the Sine-Gordon potential i.e., (9 a ) — ^/7r/8sign(/x). 
This allows us to express the dual fields in the chains as 
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Using (^), the staggered part of the transverse spin correlation functions are now given by 

x± (x,t) ee (S+ 2 S-) ~ (-lT/ a e- iBx (x 2 - vH 2 )-i (45) 

where B = ( ^1 ^ . The transverse correlations still fall off as a power-law with exponent 1/4, but they are incom- 
mensurate with a field dependent characteristic vector: 

Comparing (^5|,[l6]) with the results of Ref. ||, which corresponds to setting fi c — in the above equations, we find 
that the exponents are the same in both cases and the only effect of the chiral term is to induce a magnetic field 
dependent incommensuration. Moreover, the Pauli term does not induce any additional incommensuration in the 
staggered transverse spin correlation functions. Therefore, a varying Q should serve as a good experimental probe of 
the chiral interactions. 

In the gapless incommensurate phase of the singlet (h > h c ), the singlet excitations are described by a Luttinger 
liquid (fll]). In this case, both the gapless triplet and singlet excitations (which are probably chiral) contribute to x±- 
Here, 9 a also has a zero expectation value and the expression ( |44| ) is invalid in this phase. Instead, 

6 1>2 = -^(0s(x) ± 6 a (x)) + ^0-^- (47) 
which leads to the following asymptotic behaviour of the transverse spin-spin correlation: 



(x,t) ee (St 2 Si) ~ {-l) x/a er lB±x {x 2 -v 2 t 2 yi(x 2 - v 2 a (h)t 2 y^~ (48) 



where B± = ( P ^/^ P ) ■ Note that v is the velocity of the gapless triplet excitations and v a (h) is the velocity of the 

gapless singlet excitations. This behaviour is indeed reminiscent of a two component Luttinger liquid. The equal-time 
transverse spin correlations on the two chains fall off as a power-law with exponent K (h) = (1 + K a (h))/A, but they 
are incommensurate with characteristic vectors: 

At the critical field h — h c , since K a (h = h c ) — 2, the spin correlation function 

X ±(x,t) ~ e^ x (x 2 - v 2 t 2 )-i(x 2 - v 2 a (h c )t 2 yi (49) 



From (49) and (|45[), we see that the transition will manifest itself in the equal time spin correlation function as an 
abrupt change in the exponent i.e., it jumps from 4 to | and then varies slowly with increasing field. Such a dramatic 
change in the exponent could be easily seen in NMR and susceptibility measurements. 

For strong zig-zag interactions 7, we anticipate that the chiral term is irrelevant and that the singlet sector remains 
gapped. It remains to be seen whether the transition in the singlet space survives when fluctuations are taken into 
account. Also, these are complicated sine-Gordon theories and one should worry about solitons and kinks. An 
interesting question is whether such a physics survives in the isotropic case or when there is a finite J z coupling. 



V. CONCLUSIONS 



Chiral interactions induced by the magnetic field appear in a variety of non-bipartite quasi one-dimensional spin 
systems. We have shown that these terms could dramatically alter the physics of systems such as strongly coupled 
spin tubes and weakly coupled zig-zag ladders. The effect of these chiral interactions depend strongly on the direction 
of the applied field and could be measured experimentally. In the case of spin-tubes in a magnetic field along the 
tube, we showed that these chiral terms change the critical behaviour in the incommensurate phase, from that of a 
two component Luttinger liquid (in the absence of chiral interactions) to the usual single component Luttinger liquid. 
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However, for fields perpendicular to the axis of the tube, no chiral terms are generated and one recovers the usual two 
component Luttinger liquid. On the other hand, these terms have no effect on certain strongly coupled zig-zag and 
frustrated ladders. For the weakly coupled zig-zag ladder, we showed that the chiral term does something entirely 
novel i.e., an applied magnetic field can drive a C-IC like transition in the non-magnetic singlet sector! This results in 
gapless singlet and triplet excitations described by a two component Luttinger liquid. It is very remakable that the 
chiral term destroys the two component Luttinger liquid in spin tubes, but induces a novel two component Luttinger 
liquid behaviour in the weak zig zag ladder. These effects should be visible in measurements of various experimental 
quantities ranging from simple specific heat, NMR rates and dynamical structure factor measurements. 
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FIG. 1. Cylindrical three-leg ladder (spin-tube). 
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FIG. 2. Examples of frustrated ladders in the strong-rung coupling limit. 
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FIG. 3. Zig-zag ladder in the weak coupling limit. The arrows indicate the closed three-spins path. 
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